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Abstract  -  In  this  paper  we  describe  some  software  for  image 
reconstruction  of  electric  impedance  tomography  (EIT)  made 
on  Sharif  University  of  Technology  named  (SUT-1).  We  will 
discuss  about  image  reconstruction  methods  and  we  will  show 
the  results  of  simulation  and  real  measurement  from  SUT-1 
system. 
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I.  INTRODUCTION 

SUT-1  is  a  32  electrodes  electrical  impedance  tomograph  that 
was  designed  and  fabricated  in  Sharif  University  of 
Technology.  In  Fig.  1.  general  view  of  this  system  was 
showed.  Information  about  hardware  and  main  specification 
of  system  and  some  results  and  also  software  for  solving 
forward  problem  by  FEM  was  introduced  before  [1,2, 3, 4, 5]. 
The  main  objective  of  this  paper  is  discussion  about  some 
image  reconstruction  methods  for  EIT;  presenting  some 
results  from  image  reconstruction  software  based  on  such 
methods  and  verification  of  software  with  simulation  and 
actual  test  data  from  SUT-1. 

Fig.  1.  General  view  of  SUT-1 


II.  Methodology 


Image  in  EIT  is  mapping  of  conductivity  distribution,  by 
means  of  observation  of  electric  potential  (measurement  of 
voltage)  in  the  finite  number  of  points  in  the  surface  of  object 
that  such  potentials  made  by  injection  current  to  media.  The 
sensors  for  injection  current  and  sensing  voltages  are 
electrodes. 


1 )  Mathematical  Model : 

V.[5(P).VL(P)]  =  0  at  B  (B  is  the  object)  (1) 

5  ( P) - =  /  P  e  S  (Surface)  (2) 

dn 

\  U(P)ds  =  0  PeS  (3) 
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where  U  (P)  is  voltage  and  8(JP)  is  specific  admittance  of 
B;  in  which: 

8  (  P)  =  8  (  P)  +  y  C0£  (  P)  S  is  surface  boundary  of  B 

These  equations  are  derived  form  of  Maxwell  s  Equations 
with  some  approximation  and  an  elliptic  partial  differential 
equation  with  Neumann  boundary  condition  and  with  low 
frequency  approximation  of  current  injection  in  conductive 
media.  As  described  above  for  solving  this  equation  in 
forward  problem  FEM  is  used. 

2)  Image  Reconstruction  Methods :  Most  important  subject 
in  EIT  is  its  image  reconstruction,  because  image 
reconstruction  in  EIT  is  an  ill-posed  inverse  problem.  From 
the  time  of  advent  EIT  in  1983  very  much  works  were  made 
in  EIT  image  reconstruction  and  this  kind  of  effort  is 
continued  up  to  now.  For  reconstruction  in  mode  of  static 
imaging  a  modified  NewtonRaphson  Method  with  32 
electrodes  was  implemented.  In  this  approach  we  interested  in 
minimizing  the  function  (f>  with  respect  to  G  defined  as 
follows: 

$ =\u -vy  u -v)  (4) 

V  is  measured  voltages  and  f  is  forwad  problem  results. 

The  minimization  of  §  turns  out  to  be  the  Newton  iteration 
which  is  shown  in  this  equation,  This  is  simply  an  update  for 
g: 

Aot  =-[f'(ck)Tf\aK)Y1f\ak)T[f(aK)-V]  (5) 

For  image  reconstruction  with  this  method  some 
regularization  needed  [6].  And  also  sensitivity  method  used 
for  image  reconstruction.  It  uses  a  theorem  derived  by 
Geselowitz  [7].  The  theorem  is  best  described  using  Fig.  2  . 
When  the  conductivity  distribution  changes  from  G  to  G  +  Ag 
the  transfer  impedance  change  AZ  for  the  pairs  of  current  and 
voltage  electrodes  (A,B)  and  (C,D)  can  be  given  as: 
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r  Vw(o)  Vi’(o  +Ag) 

AZ=-  Ac? - . - c/f2  (6) 

Jn  /  / 

U  V 

where  u  is  the  potential  distribution  over  the  field  when  the 
current  Iu  is  applied  at  (A,B)  with  a  conductivity  of  a. 
Likewise  v  is  the  potential  over  the  field  when  the  current  Iv 
is  applied  at  (C,D)  with  a  conductivity  of  a  +  Aa.  This 
functional  is  used  to  help  us  solve  the  inverse  problem  in  the 
following  manner.  We  guess  a  conductivity,  (a),  and 
calculate  u  given  Iu  in  the  forward  problem.  Taking  the 
difference  between  the  calculated  u  and  the  measured  v  will 
tell  us  AZ,  using  this  we  should  be  able  to  solve  for  Aa  and 
hence  find  the  conductivity  distribution. 


Fig.  2.  Sensitivity  Theorem 


A  method  implemented  for  image  reconstructed  based  on 
Neural  Network,  another  method  is  based  on  genetic 
algorithm  direct  mapping  of  vector  of  voltage  to  conductivity 
was  tried  but  the  computation  time  for  image  reconstruction 
is  very  high,  however  some  researcher  used  these  methods  for 
image  reconstruction  in  EIT.  Other  one  is  based  on 
perturbation  method.  This  iterative  technique  was  based  on 
the  application  of  known  voltage  patterns  and  measurement 
of  current: 


This  method  was  implemented  by  Kim  [10]. 

Another  Image  reconstruction  method  is  Double  Constraint 
Method.  And  also  Layer  Stripping  Method  is  other  way.  Also 
DSM  (Direct  Sensitivity  Method)  is  suitable  for  small 
perturbation  detecting  in  conductivity.  A  very  fast  image 
reconstruction  method  with  using  one  step  of  Newton- 
Raphson  and  solving  forward  problem  by  analytical  method 
is  implemented  this  method  is  named  by  NOSER  [8].  Also  a 
fast  image  reconstruction  developed,  this  method  is  Broyden 
Quesi  Newton  Method  [9]  that  was  introduced  before.  We 
developed  a  mixed  image  reconstruction  based  on  Broyden 
and  modified  Newton  methods  that  in  first  step  of  image 
reconstruction  modified  Newton  Rophson  algorithm  is  used 
and  the  near  of  solution  Broyden  method. 


III.  RESULTS 

We  measured  voltage  and  reconstructed  images  in  the  real 
case  both  in  16  and  32  electrodes.  Fig.  3.  parts  a-f  show  some 
of  results  both  from  simulated  and  real  data.  Fig.  3a  and  3f 
illustrate  two  of  the  actual  images  using  a  simple  phantom  in 
dynamic  mode.  The  phantom  was  made  of  PVC  and 
electrodes  were  Ag/AgCl.  Star  artifact  could  be  observed  in 
the  image  resulted  from  back-projection. 

IV.  CONCLUSION 

In  this  paper,  we  discuss  about  image  reconstruction 
methods  in  EIT  and  also  some  image  reconstruction  software 
implemented  for  SUT-1  a  32  electrode  EIT  system  that  made 
based  on  such  methods  is  introduced.  Verification  of  software 
in  simulation  and  actual  data  from  SUT-1  system  was 
showed. 


For  dynamic  imaging  mode  of  operation,  image 
reconstruction  was  performed  with  16  electrodes  using  Back 
Projection  algorithm  and  iso-potential  lines.  Basically  it  was  a 
Sheffield  Algorithm  with  some  improvement  and 
modification  for  matching  with  SUT-1  system  characteristics. 
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a.  A  simulated  image  with  BP  method 
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b.  An  example  of  Broyden  Method,  Simulated 
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c.  Sensitivity  method,  simulated 


d.  Newton-Raphson  method,  simulated 


e.  The  model  of  phantom  with  p  1=150  & 
p2=  15  <Am 


f.  A  reconstructed  image  of  e  in  APT  with  16  electrodes 
using  BP(real  data) 


Fig.  3.  Some  reconstructed  images  by  simulation  and  also  real  data 


